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SUMMARY(U) 

A  number  ol  Inertial  Navigation  System  (INS)  models  ranging  from  simple 
forced  harmonic  oscillator  models  to  three-gyro,  three-accelerometer  systems 
were  studied  to  ascertain  a  compromise  between  accuracy  and  computation 
time  on  the  one  hand  and  model  complexity  suitable  for  use  in  a-vionics  system 
models  for  eff^veness  studies  on  the  other.  INS  position  and  orientation 
errors  for  various  cruise  and  acceleration  conditions  were  predicted  by  these 
models  and  the  results  shown  graphically.  The  importartce  of  various  INS 
characteristics  has  been  determined  and  the  effects  of  foe  relevant  error  sources 
have  been  isolated  and  their  propagation  in  time  plotted. 
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1.  INTRODUCTION 

The  work  described  here  is  part  of  task  AIR  86/166  to  evaluate  the  performance  of  the  Flll-C  Pave  Tack 
system.  This  task  involves  the  development  of  a  digital  computer  performance  model  which  necessarily 
encompasses  a  number  of  the  F111>C  avionics  systems,  namely  those  which  significantly  affect  navigation 
and  weapon  delivery.  Aircraft  flight  characteristics  and  the  environment  also  need  to  be  considered.  The 
model  is  inteitded  to  examine  many  possiUe  scenarios  and  options,  and  a  deterministic  approach  is  used, 
which  implies  that  raiKlom  decisions  selected  at  numerous  times  during  a  simulation  (ie  a  single  run  of  the 
performance  model)  must  be  averaged  by  repetition.  There  is  therefore  a  need  to  represent  any  system  to  be 
inoorporated  in  the  model  in  as  simple  a  way  as  possible  without  compromising  the  accuracy  of  the  overall 
trtodel.  One  such  system  is  the  Inertial  Navigation  System  (INS). 

Mathematical  models  of  generic  and  particular  INS  exist,  see  references  1  to  10.  In  simpler  models, 
references  1  to  3,  a  coiutant  drift  in  position  error  (radial)  is  cottsidered  to  be  adequate.  Overall,  models 
range  from  simple  forced  harmonic  oscillators  to  detailed  representations  of  the  gyroscopes  and 
accelerometers  that  comprise  an  INS,  incorporating  accurate  representations  of  the  sources  and  propagation 
characteristics  of  system  errors. 

Three  generic  models  ranging  from  simple  to  detailed,  were  examined  and  compared.  A  simplified 
technique  for  representing  radial  aircraft  position  error  in  an  INS  was  examined  first.  The  next  model  was 
a  two  dimensional  “tangent  plane”  system  which  predicted  the  longitudmal  and  lateral  position  errors  as 
well  as  platform  pointing  errors.  The  third  model  examined  was  much  more  complex,  incorporating  the 
effect  of  geographic  latitude  and  longitude,  altitude,  and  better  error  representation. 

The  model  is  required  to  be  suitably  accurate  while  fast  enough  to  be  used  as  part  of  the  overall  Pave  Tack 
perfonnance  model.  The  analysis  of  the  improvement  in  navigational  capability  brought  about  by  the  use 
of  Pave  Tack  requires  an  appraisal  of  the  effect  of  acceleration  on  the  INS.  This  aspect  was  critical  in 
determining  which  of  the  three  models  examined  was  most  suitable.  Software  was  written  for  the  first  two 
models  and  specific  flight  conditions  were  examined.  Comparisons  were  made  between  models  one,  two  and 
the  published  results  of  modd  3(ref.4).  It  is  proposed  to  validate  the  chosen  model  at  some  later  date, 
with  respect  to  the  effect  of  acceleration  on  INS  error,  by  monitoring  and  recording  aircraft  position  data 
obtained  from  aircraft  instrumentation  and  comparing  it  with  external  radar  tracking  data. 


2.  COORDINATE  SYSTEMS  AND  NOTATION  CONVENTIONS 

Three  right-handed  Cartesian  coordinate  systems  are  used  namely: 

(a)  True  axes  OXYZ,  are  right-handed  Cartesian  with  coordinates  of  any  point  denoted  by  (x,y,z). 
This  system  is  fixed  to  the  earth.  The  X-Y  plane  is  tangential  to  the  earth's  surface  at  the  origin,  and 
the  X-axis  points  north.  Away  from  this  base  point  the  ^tem  axes  are  not  locally  level  and  will  not  in 
general  point  north. 
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(b)  Platform  axes  OXpYpZp  fixed  to  the  aircraft,  origin  at  the  INS  origin.  The  positive  OX-axis  is 

aligned  with  the  aircraft  longitiidinal  axis  and  the  positive  Oz-axis  is  upward.  The  coordinates  of  any 
point  are  denoted  by  (Xp,yp,Zp). 

(c)  Computed  axes.  The  coordinates  of  any  point  are  (x<.,yc,2c)'  t«pre8enting  the  computed  position  or 
best  estimate  of  (x,y,z). 

The  angular  deviations  between  these  three  sets  of  coordinates  are  assumed  to  be  small.  They  can  be 
represented  simply  by  angle  vectors  as  follows: 

A  is  the  angle  vector  between  true  and  platform  axes. 

^  is  the  an^e  vector  between  true  and  computed  axes. 

Hi  is  the  angle  vector  between  computed  and  platform  axes. 

Components  of  %  are  refered  to  as  platform  errors.  Position  errors,  ie  the  differences  between  true  and 
computed  positions  are  represented  by  a  subscript  e,  (eg  x^).  Derivatives  vrith  respect  to  time  are  denoted 
by  single  or  double  dots  respectively,  over  a  symbol. 


3.  INS  CHARACTERISTICS 

In  this  section  some  phenomena  associated  with  inertial  navigation  systems  are  considered  in  order  to 
determine  their  relevance  to  the  Pave  Tack  evaluation  model. 

3.1  Schulertuning 


An  oscillatory  mechanical  system  which  is  designed  to  sense  the  direction  of  the  earth's  gravitation  is 
not  affected  by  constant  velocity.  If  however,  it  is  subjected  to  a  uniform  acceleration  it  will  in  general 
indicate  a  shift  from  the  true  direction  because  it  caiwot  distinguish  between  the  applied  acceleration 
and  that  due  to  gravity.  Consider  an  accelerometer  mounted  on  a  servo  controlled  stable  platform 
subject  to  gravitational  acceleration  g.  If  the  platform  is  tilted  by  an  angle  0  the  output  a  from  the 
accelerometer  is  integrated  twice  (to  determine  0)  and  is  fed  back  to  the  servo  to  level  the  platform.  If 
a  uniform  horizontal  acceleration  ap  is  applied  to  the  platform  the  output  of  the  accelerometer  is  given 
by 


a  =  ap  +  gsinO 


(1) 


and  the  servo  adjusts  the  platform  angle  until 


!E  +  sin0  -  ^E+e  =  0 
S  g 


(2) 
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90  that  the  pUtfonn  is  no  longer  level.  The  tilt  0  satisfies  the  relationship 


0  =  -k/  Jadtdt  =  -kgJ/Odtdt 


(3) 


where  k  is  a  constant.  The  equivalent  differential  equation  is 


S-t-kgO  =  0. 


<4) 


The  solution(ref.l6)  is  hannonic  with  period  T  and  angular  frequency  (i>  given  by 


T  =  2*/V(kg) 


and 


0)  =  V(kg). 


Since  the  damping  factor  is  zero  the  solution  is  oscillatory  and  can  be  written,  denoting  time  by  t,  as 


9  =  A  sin  (i)t. 


In  general  the  centre  of  oscillation  diffm  from  the  true  g  direction.  It  can  be  shown(ref.4,13)  that 
if  k  is  the  reciprocal  of  the  earth's  radius  R,  g,  coincides  with  g.  The  natural  period  of  oscillation, 
denoted  by  Tg,  becomes 


T,  =  2>r/(R/g) 


(5) 


and  an  applied  acceleration,  or  any  arbitrary  motion,  will  not  change  the  centre  of  oscillation  of  the 
instrument.  Taking  R  =  6371  x  10®  m  and  g  =  9.8062  m/s®  this  period  is  found  to  be  84.4  min.  The  process 
of  attempting  to  time  such  an  instrument  to  have  a  period  Tg  averaged  over  R  and  g  in  the  area  of 
operation  is  kitown  as  Schuler  tuning.  In  practice  the  equivalent  pendtilum  length  achieved  by  this 
prosess  is  usually  somewhat  less  than  R. 
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An  insight  into  this  process  is  obtained  by  noting  that  the  period  of  a  simple  pendulum  is  2xV(l/g) 
where  1  is  the  pendulum  length.  If  its  period  is  Tg  then  1  =  R.  the  radius  of  the  earth.  A  force  causing 

horizontal  acceleration  of  the  pendulum  platform  acts  on  the  point  of  support  of  the  pendulum  while  g 
acts  on  the  centre  of  gravity  of  the  bob.  If  this  coincided  with  the  earth's  centre,  as  it  would  if  1  =  R,  a 
horizontal  acceleration  of  the  platform  would  not  disturb  the  peitdulum  swing. 

A  Schuler  tuned  system  oscillates  about  true  g  with  period  Tj.  Without  Schuler  tuning,  under  a 
horizontal  platform  acceleration,  the  S3rstem  oscillates  with  period  T  about  an  erroneous  centre  of 
oscillation  g^  aiul  the  resultant  position  error  will  oscillate  with  the  same  period.  In  addition  the 
offMt  error  (resulting  from  an  accelerometer  giving  an  output  when  there  is  no  applied  lateral 
acceleratian)  is  proportional  to  t^  without  Schuler  tuning.  When  the  system  has  been  Schuler  tuned 
this  error  also  oscillates  with  period  Tg  and  is  bounded. 

3.2  The  Foucault  effect 

If  a  simple  pendulum  suspended  above  the  earth  is  started  swinging  in  a  definite  vertical  plane,  this 
plane  will  rotate  slowly  about  the  vertical  axis.  Let  Q  denote  the  angular  velocity  of  precession 
relative  to  the  earth  and  Lg  the  angle  between  the  horizontal  and  the  earth  axis  (ie  geographical 
latitude,  assuming  the  indicated  and  the  actual  gravitational  directions  are  the  same).  Then  the 
angular  velocity  of  precession  is  given  by,  see  reference  11, 


Q  =  -tugsin  Lq 


(6) 


where  tog  is  the  angular  velocity  of  the  earth. 

In  the  northern  hemisphere  the  direction  of  precession  is  clockwise  (looking  towards  the  centre  of  the 
earth)  and  anticlockwise  in  the  southern  hemisphere.  At  the  equator  Q  is  clearly  zero.  At  the  poles 
Q  =  ±<B£  and  the  pendulum's  plane  of  swing  remains  constant  in  space  while  the  earth  turns  beneath  it. 
The  period,  or  time  for  one  complete  precession  cycle,  clearly  varies  from  one  day  to  infinity.  For 
example  at  latitude  45°  the  period  will  be  '</2  days  or  34  hours.  If  the  pendulum  is  observed  with 
respect  to  a  geographic  reference  ^stem  the  projection  of  its  oscillation  on  the  x  and  y  axes  will  appear 
as  beat  waves(ref.8)  with  short  period  proportioital  to  '^(R/g)  (Schuler)  and  long  period  proportional  to 
(togsin  La)*'.  For  an  INS  system  the  Foucault  effect  on  position  errors  is  to  modulate  the  Schuler 
oscillation  at  a  frequency  O.  The  amplitude,  however,  is  second  order  and  can  be  neglected  in  this 
study.  For  platform  errors  the  Fbucault  modulation  is  a  first  order  effect(ref.6).  Fortunately  the 
platform  errors  are  of  secondary  importance  for  navigational  purposes. 

33  Geographic  latitude  and  longitude  errors 

If  the  INS  position  errors  are  resolved  into  latitude  and  longitude  rather  than  x^  and  y^,  it  can  be 
shown(ref.6)  that  latitude  error  is  of  the  form  (1  -  cos  which  is  bounded  (in  the  long  term,  eg 

over  24  hours)  and  that  longitude  error,  which  is  of  the  form. 
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[UE^  +2  (1  -cos  (i>e)'/^ 


(7) 


is  unbounded.  Howevo',  for  the  first  several  hours  the  rate  of  increase  of  the  latitude  errors  is  greater 
than  that  of  the  longitude  errors.  This  is  supported  an  analysis  of  operational  data  or  drift  rates  of 
the  AN/AJQ-20A  INS  fitted  to  the  Flll-C  aircraft(ref.I2)  which  quantifies  this  difference.  Data 
from  584  IBNS  flight  data  record  cards  were  used  for  this  purpose.  The  above  is  consistent  with  radial 
error  (Ml)  and  and  (M2)  being  unbounded  in  general,  as  foimd  in  Sections  4  and  5-  Latitude  errors 
will  also  depend  on  whether  the  direction  of  flight  is  easterly  or  westerly.  Further  consideration  is 
given  to  this  matter  in  the  next  sectioit. 

3.4  Gyroscope  aiul  accelerometer  emxs 

INS  system  errors  arise  from  lack  of  perfection  in  the  gyroscopes  and  accelerometers,  and  from  errors  in 
the  initial  values  of  INS  parameters,  such  as  alignment  direction,  position  and  velocity.  Many  errors  of 
varying  degrees  of  importance  can  be  enumerated.  They  depend  on  the  INS  configuration  (eg  the  number 
of  gyros  and  accelerometers  used)  and  the  way  in  which  the  components  are  constructed  and  combined. 
For  the  purpose  of  this  study  (following  reference  6)  gyroscope  drift  rates  ex,  Cy  and  ez  were  identified 
as  well  as  drift  rates  Cx  and  ^  due  to  mass  uifoalance.  The  latter  arises  because  of  the  inevitable  slight 
misaligrunent  of  a  gyroscope  centre  of  gravity  with  its  axis  of  support 

Of  the  many  Imown  accelerometer  errors,  bias  and  scale  factor  eTror(Te{.4)  were  identified  as  the  most 
significant  and  are  incorporated  into  M2.  Accelerometer  bias  is  the  mean  acceleration  indicated  by  the 
accelerometer  when  the  actual  input  acceleration  is  zero.  The  response  of  an  accelerometer  to 
acceleration  involves  a  scale  factor,  and  the  uncertainty  in  the  value  of  this  factor  causes  a  steady  state 
error  in  the  prediction  of  acceleration.  This  error  is  proportional  to  the  actual  acceleration,  and  is 
called  the  scale  factor  error. 


4.  FORCED  HARMONIC  OSCILLATOR  MODEL,  MODEL  M1 

The  differential  equations  derived  for  an  inertial  navigation  system  are  in  general  coupled.  However  it 
has  been  shown(ref.5),  that  if  the  coupling  terms  (eg  see  equations  (20)  and  (21))  are  neglected,  although 
the  values  of  x^(t)  and  y^ft)  obtained  are  different,  a  good  first  approximation  to  the  radial  error 
2  2 

(Xg(t)  +  yg(t))*^^  can  be  found.  The  decoupled  equations  can  be  written  in  the  form: 


Xg(t)  +  (Ogx^(t)  =  F,(t) 
ygd)  +  (agyefi)  =  Fzfi) 


(8) 
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where  tug  is  the  Schuler  frequency  and  the  Fi(t)  are  arbitrary  functions  (accelerations)  refered  to  as  forcing 
functions.  These  equations  are  seen  to  be  harmonic  osdllators.  For  each  equation  the  solution  consists  of  two 
parts,  the  comptementary  function  (Fjft)  =  0)  and  a  pardctilar  integral  dependent  on  the  forcing  function. 

These  two  parts  are  added  to  form  the  general  solution  to  the  equation.  The  comptementary  function  is 
readily  foimd  to  be 


xgft)  =  Xe(0)cos{ost  +  (l/o)s)ie(0)sin  (i>st. 


(9) 


This  is  approximated  well  by 

I 

1 


Xg(t)  =  Xe(0)  +  Xe(0)t 


(10) 


2 

since  the  Schuler  frequency  (og  (see  Section  3.1)  is  given  by  (Og  =  1.6  x  10"*  r/s,  and  therefore  tog  is  very 

small  allowing  the  approximation  0  =  sin  6.  A  number  of  cases,  distinguished  by  different  forcing  functions 
will  be  considered  and  the  results  will  be  compared  with  those  from  model  M2. 

4.1  Constant  speed 

In  this  case  the  aircraft  on  which  the  INS  is  mounted  is  flying  straight  and  level  at  a  constant  speed. 
The  forcing  function  is  assumed  to  be  a  random  variable  A(t)  with  mean  square  value  Rg  and  an 
exponentially  decaying  autocorrelation  function  of  the  form  Roe*®  1*1.  This  yields  a  particular  solution 
in  the  form  of  a  mean  square  which,  when  averaged  in  time  (denoted  by  the  vinculum),  becomes 


XRMS<*>  =  (c,at+co)’''2  (11) 

where. 


c,  =  Rq/ (<4(a^ «4^) 


and 


Cfj  =  c,(2a^/(o^  +  o^)-l/2). 
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4^  Constant  linear  acceleration 

Consider  the  application  oi  a  step  of  magnitude  a  in  acceleration  along  the  aircraft  longitudinal  axis 
for  a  duration  d  seconds  from  time  tj  to  t2  while  the  aircraft  is  flying  straight  and  level  at  constant 
speed.  The  assumption  is  made(ref.5),  that  the  forcing  function  in  the  error  is  proportional  to  the 
forcing  function  for  the  position  coordinate  x(t).  Denoting  the  unit  step  function  by  U(t),  and  the 
constant  of  proportionality  by  k,  the  forcing  functions  can  be  written: 


F,(t)  =  ka(U(t  -  t,)-U(t  -  tz)) 


Fjft)  =  0. 


(12) 


The  particular  integral  for  the  equation  in  x^ft)  is  readily  found  to  be 


x^(t)  =  — j- ((I  -  cos  <i>s(t  - 1,))  -  (1  -  cos  togCt  -  tj))!,  t>t2>t, 


(13) 


Using  trigonometrical  transformations  this  can  be  written,  using  d  in  place  of  f2  -  tj  as 


Xgft)  =  — 2"  sinfwsft- jd)Jsin(|o)sd) 
®S 


or,  since  teg  is  small. 


Xg(t)  =  kad(t-^). 


(14) 


4  J  Constant,  centripetal  acceleration 

If  while  flying  straight  and  level  at  constant  speed,  the  aircraft  commences  a  constant  centripetal 
acceleration  turn  (ie  circular  motion  in  the  horizontal  plane)  of  magnitude  ng,  the  x  and  y  components  of 
the  force  are  -ng  cos  a  and  ng  sin  a  respectively.  Where  a  is  measured  anticlockwise  from  the  y-axis. 

The  angular  velocity  is  given  by  oi  =  ng/v  so  that  the  forcing  functions  are  given  by: 
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F,(t)  =  -k,  sin(ngt/v> 


(15) 


F2(t)  =  kiCOs(ngt/v) 


(16) 


where  k]  is  the  amplitude  of  the  forcing  oscillation.  The  particular  solutions  for  this  type  of  forcing 
function  are  well  known(ref.ll),  to  be; 

Xgft)  =  cosfi*  (17) 

<0g-0i^ 


ye(t)  =  — —  sincot.  (18) 

a)g-cij2 


These  solutions  are  compared  with  those  of  model  M2  in  the  next  section.  The  constants  k,  kj  in  this 
section  and  in  Section  42  need  to  be  determined  empirically  or  by  comparison  with  a  better  model. 


5.  TANGENT  PLANE  MODEL,  MODEL  M2 

Model  M2  is  derived  in  reference  4.  There  are  a  number  of  comparable  models  in  the  literature,  references  5 
to  8.  The  error  sources  considered  in  this  model  include: 


accelerometer  bias, 
accelerometer  scale  factor, 
inherent  gyro  drift  rate, 
gyro  mass  unbalance  drift  rate, 
initial  platform  misalignment, 
initial  position  errors,  and 
initial  velocity  errors. 


In  general  the  INS  alignment  procedure  undergone  before  an  aircraft  takes  off  implies  a  correlation  between 
the  initial  position  and  initial  velocity  errors! ref .4).  This  has  not  been  included  but  will  not  affect  the 
conclusions  drawn  here. 


r 
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5.1  Potition  enor  equations 

The  radial  position  error,  (ie  the  vector  sum  of  errors  in  the  level,  nominally  horizontal,  plane)  is 
denoted  by  ig.  It  satisfies  the  equation(ret.4). 


ifg(t)  +  (i)g  i^(t)  =  -{j|£(t)  X  A(t)]r  +  [A^lr 


(19) 


where  A  w  ihe  accelerometer  signal,  Ae  is  ‘is  component  of  Ae  denoted  by  ['If.  The 

resemblance  to  equations  (8)  is  clear.  The  solutions  will  again  oscillate  with  Schuler  frequency  cog. 
The  vector  product  in  the  first  of  the  terms  on  the  right  hand  side  of  equation  (19)  is  the  error  in  sensing 
A  resulting  from  the  misorientation  The  term  Ae  includes  the  accelerometer  bias  and  scale  factor. 
This  equation  can  be  resolved  into  Xg  and  y^  components  to  obtain  the  INS  position  errors  for  a  cruising 
vehicle,  namely 


Xe(t) 


Xg(0)  cos  (ogt  +  x^(0) 


sin  ogt 
<0S 


-  Vy(0)  R(TI  -  cos  (ogt)  + 


Bx 

¥2<0)Vyt  +  —  <1  -  cos  (Ogt)  - 
“S 


r, 

“sJL 


+  EzVy  [t^  -  2/(1^  (1  -  cos  (Ogt))  -  J  Eygvt^ 


(20) 


sin  (OQt 

ye(t)  =  ye(0)  cos  (ogt  +  y^CO)  — — —  +  Vx^O)  R(1  -  cos  (ogt)  - 


-  eivx  It^  -  2/<i^  (1  -  cos  togt))  + 1  ejjgvt^ 


(21) 
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where  fi 

K,,Ky 


8 

R 


Vx.Vy 


£x#8y<*2 

Ex,Ey 


is  the  accelerometer  bias, 
are  accelerometer  scale  factors, 

is  the  acceleration  due  to  gravity, 
is  the  earth's  radius, 

are  the  components  of  v  in  the  true  reference  frame, 

are  inherent  gyro  drift  rates,  and 

are  gyro  drift  rates  due  to  mass  unbalaiKe. 


It  should  be  noted  that  dimensional  analysis  of  the  above  equations  as  presented  in  reference  4 
(equation  7.29)  indicated  discrepancies.  These  were  taken  to  be  typographical  errors.  The  equations 
shown  above  were  derived  from  earlier  work  in  reference  4  and  are  dimensionally  correct.  The  above 
expressions  can  be  evaluated  directly  to  determine  the  error  in  x  and  y  for  a  given  set  of  initial 
conditions.  The  variables  and  constants  used  in  the  evaluation  are  defined  in  Table  1  together  with 
their  initial  values.  These  equations  involve  general  acceleration  terms.  When  dealing  with  cruise 
conditiofts,  x  and  y  are  assumed  to  be  impulses  at  launch  (t  =  0)  and  zero  after  t  =  0.  To  determine  the 

position  error  when  x  and  y  are  not  constant  with  time,  the  following  differential  equations  need  to  be 
solved; 


..2  "2  2 
XgCt)  +  (iigXgd)  =  -Vy(0)g  +  Vz(0)  (y<t)  +  tOgVyt)  +  Bx  -  KxCOgVxt 


t  I  ••  2  »  I  2 

-  Eygt  +  Ezt  (y<t)  +  UgVyt]  -  EyCv  +  J  (OgVt^lg 


(22) 


ye(‘)  +  «sye(‘) 


Vx<0)g  -  Vz(0)  (x(t)  +  WgVxt)  +  By  -  KyWgVyt 


*  •  ••  2  **12 
-  Exgt  +  Ezt  [x(t)  +  tOgVxt]  +  ejjtv  +  j  a>5Vt2)g. 


(23) 


The  software  to  solve  these  equations  accoirmxxlates  general  definitions  of  x(t)  and  ^t)  in  the  form  of 
complete  sets  of  values  specified  at  the  integration  interval. 
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5^  Platform  oxar  equations 

The  INS  platform  eiTors(ref.4),  are  treated  in  a  similar  fashion.  Assuming  the  horizontal  acceleration 
v(t)  to  be  an  impulse  at  launch  and  the  velocity  v(t)  to  be  constant  for  t  >  0,  the  expressiotts  for  the 
platform  errors  are; 

Vx<‘)  =  Vx(W  +  «xt  +  Ex<v  + jwgvt^) 

Vy(t)  =  Vy{0)  +  eyt  +  E^v  + jtOgVt^) 

y^ft)  =  y^fOl  +  eit.  (24) 

In  the  cases  where  (Kt)  *  0  for  all  t  ^  0  the  following  equations  need  to  be  solved; 

..  •  2 

Vx<f)  =  ext  +  ex(v  +  ti>svt) 

Vyft)  =  eyt  +  eyfv  +  oigvt) 

yz(f>  =  «z- 

The  initial  conditions  pertaining  to  equations  (24)  and  (25)  are  also  shown  in  Table  1. 

6.  RESULTS  AND  COMPARISONS  BETWEEN  M1  AND  M2,  CRUISE 
CONDITIONS 

5.1  Variation  of  gyroscopic  drift  rate 

In  M2  the  gyroscopic  drift  rate  e'  is  assumed  to  have  a  constant  value  which  would  be  determined  by 
observation  or  from  the  INS  specifications.  It  is  shown  below  that  it  is  clearly  important  to  ascertain 
the  correct  value  of  the  gyro  drift  rate  to  be  able  to  correctly  estimate  the  errors  in  a  given  INS  system. 
M2  is  preferred  to  Ml  for  long  term  simulation  of  radial  position  error  due  to  gyro  drift  rate. 
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6.1.1  PoaiiioH  erron 

Using  model  M2,  the  variation  in  position  error  with  time  is  shown  for  four  different  values  of  the 
drift  rate,  namely  0.01, 0.1, 1.0  and  10.0  °/h,  in  figuro  1  where  logjox^Ct)  is  plotted  against  logigt.  In 
this  model  and  ty  are  assumed  to  be  identical.  The  modulation  on  these  curves  is  due  to  Schuler 

oscillation  (see  Section  3.1).  The  local  maxima/minima  in  this  modulation  occur  when  xg(t)=0. 

From  equation  (20)  it  can  be  seen  that  this  is  dependent  on  the  initial  conditions  Xg(0),  x^lO),  Vy(0), 
Vz^O)  and  the  values  of  ty,  ey.  This  accounts  for  the  different  phases  in  the  Schuler  cycle  seen 
ftom  curve  to  curve. 

No  graphs  showing  error  versus  gyro  drift  rate  were  shown  since  at  any  given  time  t]  the 
dependence  of  in  the  absence  of  this  modulation  is  linear,  ie 

Xe(£y)t=ti  =  gey/*4‘''‘^ 


where  c  is  a  constant. 

Ml  (RMS  error)  before  logarithmic  scaling  is  a 
to  the  time  axis  and  focus  at  the  point 

•] 


Taking  logarithms  to  base  10,  equation  (11)  becomes, 


The  form  of  the  corresponding  error  functions  in 
parabolic  section  (equation  (ID)  with  axis  parallel 

[CioZ-4co 
4cia  ’ 


•ogioXRMS<‘>  =  |logio<Ci“‘  +  ‘^ 


(27) 


and  since  Cg-  -jcj,  this  can  be  written 


•ogioXRMS<‘)  =  2'°gtoCi+jlogio(“*-5)^  a»>J- 


(28) 


The  graph  of  this,  using  a  log-log  scale  is  linear  in  the  time  region  of  interest.  Since  the  slope  is 
fixed  at  |,  a  reasonable  match  with  foe  results  of  M2  is  only  possible  over  the  first  several  minutes. 
An  example  is  shown  in  figure  2.  The  values  of  a  and  Rg  must  be  determined  empirically  or  by 
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error  plots  from  Ml  and  M2.  and  a  were  varied  until  a  reasonable  matdi  was  obtained.  A 
slightly  better  match  could  be  obtained  by  least  squares  fitting  but  was  oonsideied  unwarranted  for 
this  illustration. 

fi.1.2  Platform  erron 

Using  the  initial  and  operating  conditions  of  Block  7  in  Table  1  again,  the  solutions  for  platform 
errors  yx  and  yy  were  obtained.  These  are  identical  in  M2.  Error  plots  are  diown  in  figure  3  for 
both  the  cruise  condition  and  a  3g  lateral  acceleration  and  drift  rates  from  0.01  to  1.0 ’’/h.  The 
difference  between  the  accelerated  and  unaccelerated  values  were  found  to  follow  the  same  time 
curve  for  all  drift  rate  values.  The  errors  under  acceleration  x^,  and  without  acceleration  are 

related  by 


x^(t)-xe(t)  =  C(g)t 


(29) 


where  C(g)  is  a  constant  for  a  given  g.  This  effect  has  been  discussed  further  in  Section  35.  The  Ml 
model  does  not  estimate  platform  errors. 

6.2  The  effects  of  initial  velocity  error 

In  this  section  the  effect  of  the  initial  velocity  errors  ig(0),  y^lO)  on  the  position  error  predicted  by  M2 
is  examined.  Note  that  platform  error,  equations  (24)  and  (25),  is  not  dependent  on  these  initial 
velocity  errors.  In  figure  4  the  values  from  Block  1  in  Table  1  are  used  except  that  the  initial  value  of 
xgtt))  is  varied  from  0.3  m/s  to  5  m/s  and  10  m/s,  and  ei  =  ty  =  0.01  °/h  (the  ye^®'  aimlysis  is 
identical).  The  three  position  error  plots  are  shown  agaiitst  time  on  a  log-log  scale.  For  the  first  40  min 
there  is  no  correspondence,  but  after  that  the  curves  approach  each  other  asymptotically. 

65  Initial  misaligiunent  error 

Initial  misaligrunent  error  terms  yx(0)/  yy(0),  yz(0>  M2  propagate  into  both  the  position  and 
platform  error  predictions. 

6  J.  1  PooitioH  error 

Three  values  of  misalignnrent  error  10, 100  and  1000  arcsec  (ie  0.05,  05  and  50.0  irvad)  are  used 
together  with  a  gyro  drift  rate  of  0.01  °/h.  (See  Table  1,  Blocks  1,  3  aird  4  for  the  parameter  values 
used).  The  effect  of  an  initial  error  in  aligiuttent  is  seen  to  be  significant  and  long  Iasting,see 
figure  5. 
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63  J  Platform  error 

It  can  be  seen  from  the  platform  error  equations  (24)  that  the  initial  alignment  errors  are  simply 
additive  constants  and  effect  the  origin  but  not  the  ^pe  of  the  linearly  increa^g  platform  errors. 
The  values  of  1(X)  and  1000  arcsec  are  very  large  compared  to  the  irrcrease  in  platform  error  with 
tirtre  (eg  50  s).  Since  graphs  of  these  quantities  would  eff^vely  be  strai^t  lines  they  have  not 
been  plotted.  The  10  arcsec  case  is  discussed  uruler  Section  7.1.2  and  a  graph  of  the  corresponding 
platform  error  versus  tin«  is  shown  for  refererrce  with  platform  errors  under  acceleration. 


7.  RESULTS  AND  COMPARISONS  M1-M2,  UNDER  ACCELERATION 

Both  longitudinal  (along  aircraft  roll  axis)  lateral  (along  aircraft  pitch  axis)  and  centripetal  accelerations 
were  examined.  In  Ml  and  M2,  lateral  and  longitudinal  acceleration  give  rise  to  errors  of  the  same 
magnitude. 


7.1  Constant  linear  acceleration 

7.1.1  Poeition  error 

The  effects  of  constant  acceleration  (longitudinal  or  lateral)  varying  from  03  g  to  3  g  are  shown  in 
figure  6.  The  error  without  acceleration  and  the  differences  between  the  two  are  shown.  They  are 
oiily  plotted  for  50  s  as  this  is  long  enough  to  indicate  the  effects  of  such  accelerations. 

A  comparison  of  radial  errors  for  Ml  and  M2  under  constant  acceleration  is  shown  in  figure  7. 
Although  the  error  in  M2  is  clearly  non-linear  the  Ml  error  is  given  by  equation  (14)  which  is 
linear.  The  value  of  k  for  this  model  has  to  be  determined.  The  least  squares  technique  was  used  to 
ascertain  its  value  for  figure  7  where  an  acceleration  of  3  g  was  applied  for  50  s  (k  =  4.49  x  lO*^). 

7.13  Platform  error* 

In  figure  8  platform  errors  are  shown  for  linear  accelerations  of  03  g  to  3  g.  Drift  rates  are  seen  to  be 
from  2D  x  KT*  to  3.1  x  Iff^  r/s  or  24  to  38  arcsec/min  compared  to  213  atcsec/min  for  the  postulated 
unaccelerated  INS. 

73  Constant  centripetal  acceleration 

73.1  PotitioH  error* 

A  series  of  graphs  illustrating  the  effect  of  constant  centripetal  acceleration  are  ^own  in  figure  9. 
The  reference  graph  for  the  unaccelerated  case  is  also  shown,  as  are  the  differences  between  the 
accelerated  and  unaccelerated  cases.  Oscillation  of  the  error  as  the  platform  is  flown  through 
successive  complete  circles  can  be  seen.  Using  the  relationship  for  radial  velocity  v^/r  >  ng,  where  r 
is  the  radius  of  turn  and  n  is  the  lateral  acceleration  (constant)  in  multiples  of  the  acceleration  due 
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to  gravity,  g,  the  turn  duration  is  given  by  t  -  s/v  (where  s  is  distance  around  the  circular  path)  and 
the  duration  of  one  complete  turn  is  readily  shown  to  be  t  =  2Kv/ng.  For  the  3  g,  v  =  300  m/s  case 
used  for  figure  9,  t  ^  64.1  s.  The  error  is  seen  to  oscillate  with  fids  period. 

Equations  (17)  and  (18)  representing  position  error  in  Ml  under  centripetal  acceleration  are  again 
linear  over  the  time  interval  shown  in  figure  9.  A  graph  of  this  is  not  shown.  The  slope  could  be 
determined  empirically  or  by  comparison  with  a  better  model,  eg  M2.  The  maximum  difference  in 
predicted  error  between  Ml  and  M2  would  be  half  fiie  maximum  error  between  the  M2  graph  and 
the  reference  (unaccelerated)  case. 

7.2.2  Platform  errort 

It  can  be  seen  from  equations  (25)  that  the  platform  errors,  dependent  only  on  v  are  the  same  for 
both  the  linear  and  the  centripetal  acceleration  cases  (figure  8). 


8.  COMPARISON  OF  MODEL  M2  WITH  MORE  DETAILED  MODELS 

Possible  improvements  to  M2  include  the  use  of  a  reference  axes  system  based  on  latitude  Lq,  longitude  Lg 
and  altitude  A.  A  number  of  models  which  use  the  OL^LoA  frame  of  reference  were  considered(ref.4  to  8) 
and  an  attempt  was  made  to  consider  the  trade-off  between  more  accurate  repuesentation  of  the  INS  errors 
and  increased  complexity /computing  time  of  these  models  compared  to  M2.  In  reference  6  the  expression 
derived  for  the  error  in  latitude  in  re^nse  to  white  noise  is 


(it- jsin2Xt) 
73? 


(30) 


where  N  is  a  constant  gyro  drift  rate  power  spectral  density  function  (••  1  arc  min^/h^)  and  i  is  celestial 
longitude  rate.  Figure  10  shows  the  RMS  error  (the  square  root  of  the  expression  in  equation  (30))  as  a 
function  of  time  for : 

( 1 )  westerly  flight  at  latitude  35°  and  300  m/s, 

(2)  at  earth  rate  (15  °/h),  the  stationary  case  (Lo= 0 ),  aitd 

(3)  easterly  flight  at  latitude  35°  and  3(X)  m/s. 

The  difference  between  these  cases  starts  to  become  significant  after  about  two  hours,  indicating  that  a 
model  that  does  not  take  the  geographical  direction  into  account  should  not  be  relied  on  after  this  length  of 
time. 
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When  the  INS  platform  is  su^ect  to  vertical  acceleration  and/or  vertical  velocity,  the  computations  are 
inherently  unstaUe  because  of  the  variation  of  gravity  with  altitude(ref.7).  A  wei^ted  combination  of 
altimeter  and  INS  computed  altitude  can  be  used  to  deal  with  this  problem.  Models  incorporating  this 
effect  have  been  examined(ref.4  to  8).  They  were  considered  to  be  beyond  the  scope  of  the  currott  task  for 
which  the  INS  model  was  intended,  namely  to  assist  in  the  evaluation  of  the  Pave  Tack  system. 

These  considerations  iivdicate  that  M2  wotild  become  inadequate  after  2  hours  of  ^ulated  run  time,  and 
would  have  to  assume  a  knowledge  of  altitude  from  some  other  source.  From  this  examination  of  OLqLqA 
based  models  it  was  considered  that  the  increase  in  computational  complexity  and  the  corresponding 
increase  in  computing  time  compared  to  M2  were  not  warranted. 


9.  CONCLUSIONS  AND  RECOMMENDATIONS 

The  characteristics  of  INS  systems  including  gyro  and  accelerometer  errors,  initial  alignment  and  other 
initialisation  errors,  the  Schuler  and  Foucault  cycles  and  the  altitude  prediction  problem  have  been 
considered  for  inclusion  in  an  INS  nwdel  stiitable  for  use  in  the  evaluation  of  the  F111-C  Pave  Tack  system. 
Three  types  of  model  varying  from  simple  forced  harmonic  oscillators  to  detailed  three-gyro, 
three-accelerometer  systems  were  considered  and  compared. 

A  “tangential  plane”  model  M2  was  found  to  give  suitably  reliable  predictions  of  INS  position  and 
orientation  errors  for  up  to  two  hours  flight  time  under  cruise  conditions,  while  satisfying  time  and 
complexity  restrictions  imposed  by  its  intended  role  in  an  Ftll-C  Pave  Tack  effectiveness  model.  In 
addition  the  model  was  able  to  predict  the  effects  of  linear  and  centripetal  accelerations  on  these  errors. 
INS  predictions  for  a  range  of  cruise  aivl  acceleration  conditions  have  been  produced  for  use  in  the  Pave  Tack 
task.  The  rates  of  propagation  of  various  types  of  error  with  time  were  also  demonstrated  graphically.  No 
empirical  data  were  available  to  check  the  predictions  of  errors  under  acceleration,  nor  were  they 
generally  available  from  any  published  theoretical  analyses  or  manufacturers'  specihcations. 

It  is  recommertded  that  suitable  records  of  INS  output  and  aircraft  flight  position  from  external  tracking 
instruntents  be  obtained  during  Flll-C  trials  to  validate  the  model,  particularly  the  errors  under 
acceleration. 
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I  TABLE  1.  CRUISE  AND  ACCELERATION  PLATFORM  PARAMETERS 

I  The  teim  Block  is  used  to  refer  to  sets  of  input  data. 

( 

The  following  parameters  ate  used  for  all  Blocks:  Cx  runs  through  0.01,  0.1, 1.0,  10.0,  Cy  =  Cx,  Cz  =  0.01, 
speed  =  400  m/s. 

For  Blocks  1  to  5,  acceleration  is  0  m/s^. 


Initial  Velocity  Error 

XgfO),  ye(0) 

Misalignment  Error 

VxW), 

03  m/s 

5  m/s 

10  m/s 

10  arcsec  (0.05  mrad) 

BLOCK  1 

BLOCaC2 

BLOCK  3 

100  arcsec  (5  mrad) 

BLOCK  4 

1000  arcsec  (50  mrad) 

BLOCK  5 

'  BLOCK  6 

I  ACCELERATION  RUNS  THROUGH  05  g,  1  g,  2  g,  3  g. 

ejc  =  ejf  =  ei  =  0.01  ®/hr,  s^*ed  =  400  m/s 

ie(0)  =  ye(0)  =  03m/s 

yx(0) »  yyfO) »  O.OS  mrad _ 

<  BLOCK  7 

Acceleration  3  g  (lateral  or  longitudinal  as  required) 

Cx,  Cy  run  through  0.01, 0.1, 1.0, 10.0,  ez  =  0/)1 

XgfO)  =  YgiO)  =  03  m/s 

Vx<0)  ^  Vv(O)  =  0.05  mrad _ 
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A 

Bx.By 

C,Co<Cl 

k 

Kx/Ky 

I 

R 

g 

t 

V 

Vx-Vy-Vj 

x,y,z 

*c'yc'^c 

xe'Xe'^e 

2a 

EX/  Ey»  Ez 

EX/  Ey  ,  Ez 

e 

♦x'  ♦z 
yx-Vy.Vz 


NOTATION 

Acceleration  signal  vector 
Components  of  accelerometer  bias 
Constants 
Constant 

Components  of  scale  factor 
Position  error  vector 
Earth  radius 

Acceleration  due  to  gravity 
Time 

System  velocity 

System  velocity  components,  true  coordinates 

True  position  coordinates 

Computed  position  coordinates 

Errors  between  true  and  position  coordinates 

The  angle  vector  between  true  and  computed  axes 

Constant  gyro  drift  rate  components 

Gyro  drift  rate  components  due  to  mass  unbalance 

general  angle 

celestial  longitude 

Components  of  the  angle  vector  between  true  attd  platform  axes 
Components  of  the  angle  vector  between  computer  and  platform  axes 
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»>o-v 

“E 

(0$ 


(Of.  The  angular  rates  of  true,  platform  and  computed  axes  with  re^Tcct  to  inertial  space, 
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Figure  3.  The  variation  of  platform  error  with  gyro  drift  rate 
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Figure  6.  The  effect  of  linear  acceleration  on  position  error 
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